Measurement of gravitational spin-orbit coupling in a binary pulsar system 
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In relativistic gravity, a spinning pulsar will precess as it orbits a compact companion star. We 
have measured the effect of such precession on the average shape and polarization of the radiation 
from PSR B1534+12. We have also detected, with limited precision, special-relativistic aberration of 
the revolving pulsar beam due to orbital motion. Our observations fix the system geometry, including 
the misalignment between the spin and orbital angular momenta, and yield a measurement of the 
precession timescale consistent with the predictions of General Relativity. 
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1. INTRODUCTION 

Parallel transport of the angular momentum of a gyro- 
scope moving in curved spacetime leads to geodetic pre- 
cession pj . In the solar system, the only observed exam- 
ple is precession of the Earth-Moon system as it orbits 
the Sun 0- The recently launched Gravity Probe B ex- 
periment plans to measure the geodetic precession of a 
gyroscope in Earth orbit as well as its gravitomagnetic 
Lense-Thirring precession Q . 

Binary pulsar systems are important laboratories for 
gravitational physics , in part because the strong self- 
gravity of neutron stars (GM/Rc 2 ~ 0.2) raises the pos- 
sibility of deviations as large as order unity in some al- 
ternate gravity theories — even theories that agree with 
general relativity (GR) in weak-field tests ||. Imme- 
diately after the discovery of the first binary pulsar, 
PSR B1913+16 6], it was realized that geodetic pre- 
cession could lead to variations in the path of the ob- 
server's line of sight across the pulsar's magnetic pole, 
and hence changes in the radiation pattern at Earth 0. 
Recently, the anticipated pulse-profile variations have 
been observed, in good qualitative agreement with pre- 
dictions 0, 0, , but uncertainties in the intrinsic beam 
shape prevent a quantitative measurement of the preces- 
sion rate. 

Here we describe observations of the doublc-ncutron- 
star system PSR B1534+12, in which precession is chang- 
ing the observed pulsar profile by about 1% per year 
[ill fl2j|. We outline a new, general technique for com- 
bining measurements of long-term quasisecular varia- 
tions caused by geodetic precession with those of pe- 
riodic variations induced by special relativistic aberra- 
tion modulated by orbital motion. Together, the obser- 
vations allow a quantitative estimate of the precession 
rate, independent of the unknown pulsar beam shape. 
Although the precision is still low, we show that this 
model-independent precession rate is consistent with the 
predicted rate in GR. The polarization properties of the 




i 

Orbital plane | Magnetic pole 

V 

To observer 



FIG. 1: The spin and orbital geometry we derive for 
PSR B1534+12, after Q. The vectors I and J denote the 
plane of the sky, while i and j show the plane of the pulsar 
orbit. The pulsar spin axis at the current epoch is S and the 
orbital angular momentum direction is k; the precession cone 
with opening (misalignment) angle 5 is shown. The angle be- 
tween S and the line of sight is the supplemental angle 
A = n — £ is shown. The angle between S and the magnetic 
pole is a, and (5 — £ — a is the minimum impact angle of 
the magnetic pole on the line of sight. A second cone indi- 
cates the sweep of the magnetic pole at the current epoch. 
The projection of S on the plane of the sky (indicated by 
the dashed line) provides angle rj, measured counterclockwise 
from the ascending node; —n is shown. Spherical geometry 
gives cos 5 = — sin i sin A sin rj + cos A cos i. 



pulsar signal are also changing: with the additional as- 
sumption that the magnetic field structure is dipolar, we 
can determine the angles between the spin and orbital 
angular momenta and the line of sight to the pulsar. 

The spin and orbital geometry is shown in Fig. The 
pulsar spin axis S will precess around the total angu- 
lar momentum, which is very well approximated by the 
orbital angular momentum direction k. This precession 
will cause potentially detectable periodic variations of the 
projection of S on the plane of the sky, and on the inch- 
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nation of S with respect to the observer. In GR, the time 
averaged precession rate of the pulsar can be written : 
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where G is Newton's constant, M Q is the mass of the Sun, 
c is the speed of light, mi and m 2 are the pulsar and com- 
panion masses, respectively, Pb is the orbital period, and 
e is the eccentricity. A description in generalized theories 
of gravity is given in Q. For PSR B1534+12, using the 
stellar masses determined through high-precision timing 
[l3|. the precession rate predicted by GR is 0.51° /yr. 



2. OBSERVATIONS AND ANALYSIS 

Observations were made with the 300-m Arecibo radio 
telescope, using the "Mark IV" data acquisition system 
[lij at an observing frequency of 430 MHz. The signal 
was processed using coherent dedispersion, providing full 
polarization information as well as a pulse shape unaf- 
fected by dispersive smearing in the interstellar medium. 
Data acquisition details have been described elsewhere 

The data span the interval from mid-1998 to mid-2003, 
incorporating some 400 hours of observing time. The 
pulsar was observed biweekly or monthly. The observed 
signal strength varied widely because of interstellar scin- 
tillation; only epochs with high signal-to-noise ratio were 
used here. Campaigns of roughly 12 contiguous observ- 
ing days were also conducted every summer except 2002. 

The cumulative pulse profile from the representative 
2001 June epoch is shown in Fig. [21 Here we also show 
a fit of the linear polarization to the standard "rotating 
vector model" (RVM) [lj|, in which the position angle 
of linear polarization ip is assumed to be parallel to the 
plane of curvature of magnetic dipole field lines rotating 
with the star, giving 



tan [tp O) - i'o] = 



sin a sin (<f> — 0q) 



cos a sin £ — sin a cos £ cos (cf> — 4>o) ' 

(2) 

where 4> is the pulse phase, 4>q and ipo are constants, a is 
the magnetic inclination angle, and ( is the angle between 
S and the line of sight. 

The position angle sweep is observed over most of the 
pulsar period, and the model fit is generally good. There 
are strong deviations from the model near the pulse peak, 
as is often seen for "core" profile components ^(|. We 
exclude this region from our fits. The data are con- 
sistent with a roughly orthogonal rotator model, with 
a = 102.8 ± 0.5° and the line of sight passing between 
the magnetic pole and the stellar equator, within a few 
degrees of the magnetic pole. We also consider the time 
evolution of a and the impact parameter of the line of 
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FIG. 2: Top panel: the position angle of linear polarization 
in 2001 June, measured clockwise on the plane of the sky (the 
convention in Q), with best fit rotating vector model (RVM) 
overlaid. Only the position angle points indicated by large 
dots were used in the RVM fit; these were weighted by their 
uncertainties, with a small uncertainty added in quadrature 
to account for deviations from the RVM. Middle panel: total 
intensity (solid) and linear polarization (dashed) profiles in 
2001 June. This profile is very similar in shape to our "ref- 
erence" profile Pq. Inset: evolution of impact angle f3 with 
time. Bottom panel: "Difference" profile Pi, representing es- 
sentially the time-derivative of the observed profile. 



sight on the magnetic pole f3, using cumulative profiles 
from each campaign, one especially strong biweekly ob- 
servation, and an earlier coherently-dedispersed profile 
from observations with the "Mark III" data acquisition 
system 0. As expected, the data are consistent with no 
evolution of a. However, as shown in Fig. [21 (3 is changing 
with time, at a rate df3/dt = —0.21 ± 0.03°/yr, resulting 
in a larger impact parameter at later times. This change 
in (3 is direct evidence of geodetic precession, and can be 
related to the system geometry and precession rate in a 
simple fashion df3/dt = f^ pin cos 77 sin i (see Fig.^l. 

The shape of the profile is also changing with time 
and orbital phase, allowing a completely independent 
probe of the precession. Secular changes in the profile 
were first noticed at 1400 MHz [ll|, but evolution of the 
430-MHz emission only became apparent with coherently 
dedispersed observations 0] . Shape variations are more 
difficult to connect directly to the precession rate than 
polarization changes. It can be done with an assumed 
model of the beam sha pe, as has been attempted for 
PSR B1913+16 0, ES [3, but in that case the true 
beam shape is still debated and the results are therefore 
less than satisfactory. 

Here we note that it is possible to make a model- 



3 



independent precession estimate by also measuring the 
orbital modulation of the profile shape caused by 
aberration — a special-relativistic effect independent of 
strong-field gravity. Aberration shifts the observed an- 
le between the line of sight and spin axis by an amount 



5a( = — - r [" cos?? S (u) + cos i sin 77 C (u)] , (3) 
smi 



where f3\ = nx/y/[l — e 2 ) is the characteristic velocity 
of the pulsar, with the orbital frequency n = 2ir/Pb, the 
projected semimajor axis x = aisini/c, and the eccen- 
tricity e all available from timing data, and where C (it) = 
cos[uj+ A e (u)]+e cos lu and S(u) = sin[uj+A e (u)}+esmui 
are functions of the time-dependent angle of periastron 
passage u) and the eccentric anomaly u through the true 
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anomaly A e (u) = 2 arctan 
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Now let P(C) be any function defined on the observed 
pulsar signal that depends on the viewing angle (such as 
integrated intensity, component width, polarization frac- 
tion, etc.). For small changes in the impact parameter, 
we Taylor-expand P(C) as P(Co) + CP' where prime de- 
notes derivative with respect to (. The effects of aberra- 
tion and precession can then be written 
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The unknown beam shape enters only through F', which 
can be eliminated by dividing these two equations. Mea- 
surements of both the orbital variation of F and its sec- 
ular drift thus allow tan r\ and S7^ pm to be determined in 
a model independent way. 

For PSR B1534+12, we have measured the evolution 
of the total intensity profile (Fig. ■ The strongest data 
scans from the annual campaigns were averaged into 12 
orbital phase bins, and analyzed together with the strong 
biweekly scans. Profiles with unusually low signal-to- 
noise ratios or suspect calibrations were discarded. We 
used standard principal component (PC) analysis tech- 
niques 0] to derive orthogonal "reference" (Po) and "dif- 
ference" (Pi) profiles that completely described the pro- 
file evolution; Pi and a single-epoch profile very simi- 
lar to P are shown in Fig. The pulse profile P of 
each observation is well modeled as a linear combination 
P = coPq + ciPl. Because the overall amplitude var- 
ied with scintillation, we chose as our observable quan- 
tity F the ratio c\/cq. The PC analysis provided es- 
timates of Co, ci, and their uncertainties, which we in- 
dependently checked through a frequency-domain cross- 
correlation technique in which a linear combination of 
the two profiles was fit in an iterative manner. Finally, 
we simulated the cross-correlation analysis to assess its 
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FIG. 3: The shape parameter F — ci/co (see text) is shown 
as a function of date in the main panel and aberration phase 
(essentially the true anomaly corrected for the advance of pe- 
riastron) in the inset. The best-fit model is shown by the 
solid line in each panel, and in the orbital- phase plot, the GR 
prediction based on the RVM model is indicated by the dot- 
ted line. We have included a small error in quadrature with 
the measurement errors to account for smearing caused by 
averaging over a range of orbital phase. To account for sys- 
tematic errors, variable data quality, and uneven time sam- 
pling, we have used a bootstrap analysis [ijj to estimate the 
uncertainties on model parameters. The resulting values and 
uncertainties are in good agreement with estimates obtained 
by scaling measurement errors to obtain a reduced-^; 2 of 1. 



sensitivity to systematic errors induced by imperfect cal- 
ibration or polarization cross-coupling, finding that such 
problems should be negligible in our dataset. 



The secular trend in F is evident in Fig. |21 This cor- 
responds to a decrease in the intensity of the core region 
of the profile relative to the lower level emission in the 
wings, consistent with early indications from 1400 MHz 
data, and as expected if precession is moving our line 
of sight away from the magnetic pole, as indicated by 
the polarization analysis above. The residuals after re- 
moving the best fit line are shown as a function of orbital 
phase. A simultaneous linear fit of F as a function of date 
and of S(u) and C{u) gives the constraints Qf m sin 2 i = 



0.421^5 /yr and flf™ sin i tan i cot 77 = 0.421°;^ /yr 
(68% confidence), where we have used (3\ — 0.67 x 10~ 3 . 
A fit of F only as a function of date yields a x 2 value 
that is 15% higher than that for the full fit. 
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3. DISCUSSION 

We have observed both long- and short-term varia- 
tions in the pulse shape of PSR B1534+12, as expected 
from geodetic precession and aberration. Assuming a 
dipolar field geometry and GR, the impact parameter 
change d/3/dt = flf™ cos 77 sin i = -0.21 ± 0.03°/yr 
yields a measurement of the previously unknown an- 
gle T] = ±115.0 ± 3.8°. Using the pulse timing value 
sin 7 = 0.975 | l3| but making no assumptions about the 
validity of the RVM, our measured pulse profile varia- 
tions yield, using equations 4 and 5, the consistent re- 
sult 77 = ±103 ± 10°. Moreover, we may now solve 
for the precession rate, n\ p ' m = 0.44±g;^7yr (68% 

confidence) or flf™ = 0.44±^ 4 °/yr (95% confidence). 
This value compares well with the GR-predicted rate of 
^ pin = 0.51°/yr. 

The misalignment angle S between the spin and or- 
bital angular momenta can also be constrained. The 
angle A is known from the polarization studies. Only 
the absolute values of sin 77 and cos i are known, but our 
profile fit requires that cos i tan 77 > 0. Therefore there 
are two possible geometries: i — 77.2° and 77 = —115°, 
which gives 8 = 25.0 ± 3.8°, or i = 102.8° and 77 = 115°, 
which gives S = 155.0 ± 3.8°. Both give identical results 
for precession in GR. As the angular momenta were al- 
most certainly aligned before the second supernova, the 
smaller misalignment value is favored on astrophysical 
grounds |2(i| . The preferred geometry then has i = 77.2°, 
77 = -115.0 ±3.8°, and 5 = 25.0 ±3.8° (FigHJ. The mis- 
alignment angle can be used to constrain mass loss and 
asymmetry in the second supernova. A full analysis will 
be published elsewhere, along with a study of the two- 
dimensional beam geometry of PSR B1534±12. 

Although the precision is, as yet, limited, this is the 
first beam-model-independent measurement of the pre- 
cession rate of a binary pulsar and reconstruction of the 
full three dimensional geometry of a binary pulsar sys- 
tem, including the misalignment angle. Future prospects 
for improvement include direct estimation of the angle 77 
by combining scintillation studies |2lj | with polarimetry 
y. The derived geometry also allows us to predict the 
effects of aberration on the pulse timing |4|; this will in 
principle allow more precise timing tests of GR in future. 

We emphasize that the general technique of combining 
observations on the orbital and precessional timescales 
to make model independent precession rate estimates 
is potentially far more general than the particular ex- 
ample given here. An especially interesting prospect 
is the recently discovered highly relativistic system 
PSR J0737-3039 H3; with ft nearly twice as large as 
that of B 1534+12, and a predicted precession timescale 
for the recycled pulsar of only 75 years, both effects will 
be quickly measured for this new system. 
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